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ABSTRACT 

In the low-redshift Universe, the most powerful radio sources are often associated with 
gas-rich galaxy mergers or interactions. We here present evidence for an advanced, gas- 
rich (‘wet’) merger associated with a powerful radio galaxy at a redshift of z^2. This 
radio galaxy, MRC 0152-209, is the most infrared-luminous high-redshift radio galaxy 
known in the southern hemisphere. Using the Australia Telescope Compact Array, we 
obtained high-resolution CO(1-0) data of cold molecular gas, which we complement 
with HST/WFPC2 imaging and WHT long-slit spectroscopy. We hnd that, while 
roughly Mh2 ~ 2 X 10^° Mq of molecular gas coincides with the central host galaxy, 
another Mh 2 ~ 3 x 10^° Mq is spread across a total extent of ~60 kpc. Most of this 
widespread CO(l-O) appears to follow prominent tidal features visible in the rest- 
frame near-UV HST/WFPC2 imaging. Lya emission shows an excess over Hell, but 
a deficiency over Lir, which is likely the result of photo-ionisation by enhanced but 
very obscured star formation that was triggered by the merger. In terms of feedback, 
the radio source is aligned with widespread CO(1-0) emission, which suggests that 
there is a physical link between the propagating radio jets and the presence of cold 
molecular gas on scales of the galaxy’s halo. Its optical appearance, combined with 
the transformational stage at which we witness the evolution of MRC 0152-209, leads 
us to adopt the name ‘Dragonfly Galaxy’. 

Key words: galaxies: individual: MRC 0152-209 - galaxies: high-redshift - galaxies: 
active - galaxies: interactions - galaxies: starburst - ISM: jets and outflows 


1 INTRODUCTION 

A crucial epoch in the star formation history of the Uni¬ 
verse occurs at 1 < z < 3, where the evolution of the star 
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formation rate density peaks and the assembly of present- 
day galaxies is in full progress (e.g. Madau et al.jlQQ^ Lilly 


et al.|19% Pentericci et al.|2001| [Hopkins|2004| [Hopkins fc 


Beacom 2006 Ravindranath et al.|2006 Lotz et al.|2006[ see 


Shapley 2011 for a review). This is also the era in which the 
quasar density shows a maximum (e.g. Schmidt et al. |I^ 
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Richards et al.|2006| Aird et al.|2010 1, indicating that there 
is a strong link between the evolution of massive galaxies 
and the growth of their central super-massive black holes. 

While numerical simulations have addressed important 
aspects of the evolution of massive high-z galaxies, such as 


hierarchical merging and gas accretion (e.g. 

Springel et al. 

2006 

Dekel et al.||2009| Bournaud et al.||201] 

), historically. 


our observational knowledge about galaxy assembly pro¬ 
cesses at z > 2 has been dominated by optical and IR studies 


(e.g. Shapley 20111. While such studies provide important 


information on the star-formation history of galaxies, they 
often lack insight into the source that will drive future star 
formation. For this, we need to investigate the fuel for star- 
burst and AGN activity, namely cold molecular gas. 

To date, molecular gas has been detected in a large 
number of galaxies at z > 1, mostly through unre solved de- 
tections of the strongest tracer ^^CO(J, J-1) (see Carilli & 
Walter|2013 Genzel et al.|201^ and references therein). Sig¬ 


nificant progress in our understanding of molecular gas at 
high redshifts has come with the introduction of wide-band 
receivers at facilities that operate below 50 GHz, such as the 
Australia Telescope Compact Array (ATCA) and Karl G. 
Jansky Very Large Array (VLA). These facilities make it 
possible to observe the lowest-J transitions of CO, which 
are among the most robust tracers of the overall molecular 
gas content, including any low-density and subthermally- 


excited components (e.g. 

Papadopoulos et aL|2000l|Hainline 

et al.|l2006l 

Carilli et al. 

|2010| jlvison et al.||2011| Riechers 

et al.l|2011a 

Swinbank et al. [2011 1. By studying the CO- 


excitation down to the lowest J-transitions, indirect results 
were obtained regarding the presence/absence of widespread 
molecular gas reservoirs in high-z galaxies (e.g. Daddi et al. 


2008| [Pannerbauer et al.||2009| pVravena et al.|[^10| [Riech- 

ers et al. 2011b| Danielson et al. 201 1| Daddi et al.||2014 |. 

Despite this, direct results on the distribution of molecu¬ 
lar gas in intermediate- and high-z galaxies are still rel¬ 
atively sparse, although they are growing in number. For 
example, regularly rotating CO discs have been imaged on 
scales of up to ~ 14kpc in various types of high-z galaxies 


(Genzel et al.||2003| 

2013 

Daddi et al.||2010| Tacconi et al. 

20101 

2013 podge et al. 

2012 

), while gaseous tidal debris 


et al.|2008 Engel et al.|2010| an d two QSOs | |Riechers et al. 
2008 Villar-Marti'n et al.||2013a K 


Among the best studied high-z objects are high-redshift 
radio galaxies (HzRGs, with Lsoomhz > 10^^ WHz e.g. 


Rbttgering et al. jl9^ . HzRGs are among the most massive 


and active galaxies in formation, and often reside in galaxy 
overdensities, as expected from progenitors of current-day 


giant central-cluster ellipticals (e.g. Rocca-Volmerange et al. 


2004| jVenemans et al. [20071 Seymour et al. 

2007 Miley & 

De Breuck|[2008 De Breuck et al. 2010 Fal 

er et al.|[2010 

Barthel et al. 2012| Wylezalek et al. 2013 

Drouart et al. 

2014 Hatch et al.[ 2014[ [Dannerbauer et a 

[2014[ Podiga- 


terns where molecular gas is spread across 30-40 kpc scales 


(|Emonts et al 

2013 

), across multiple components (e.g. De 

Breuck et al. 

20^ 

Ivison et al. 2008), in IR-bright com- 

panion galaxies ( 

Ivison et al.[[2012) or in the galaxy’s halo 

(Nesvadba et al. 

2009 Emonts et al. 2014). In addition, a 


(20041 and Emonts et al. (20141 through alignments seen 


between the CO and the radio jet axis (Sect. |4d^ , which re¬ 
semble jet-alignments seen in optical/UV, X-ray and submm 


emission (McCarthy et al.jl987 Chambers et al 

11987 

Car- 

|illi et al.|2002 

Stevens et al.[[2003[ 

Small & Blundell||20131. 

Considering these CO results. 

Ivison et al. 

(2012 1 pro- 


posed that gas-rich galaxy mergers are ubiquitous among 
starbursting HzRGs. However, the existing CO studies suf¬ 
fer from a low spatial resolution (often many tens of kpc), 
hence accurate knowledge about the properties and distri¬ 
bution of molecular gas in HzRGs remains crucially lacking. 
It is also not clear to what extent other processes, such as fil¬ 
amentary cold-mode accretion of gas from the environment 
( [Keres et al.||2009 Dekel et al.||2009j Ceverino et al.pOlO l, 

may contribute to the gas supply and triggering of activ¬ 
ity. Therefore, detailed imaging of the cold molecular gas 
in HzRGs is essential to further investigate to what extent 
their starburst and AGN activity is linked to gas-rich galaxy 
mergers. 

In this paper, we present a study of radio ga laxy MRC 
0152-209 at z = 1.92 (Pi. 4 GHz ~ I x 10^® W/Hz; |Pentericci 


et al.|[2000 |. MRC 0152-209 is the most IR-luminous HzRG 


known in the southern hemisphere (Drouart et al. 2014). 


From a detailed analysis of its IR spectral energy distribu¬ 
tion, Drouart et al. (2014) found that MRC 0152-209 has a 
starburst-IR luminosity of Lpm- SB ~ 1.8 X lO^^L©, which 
is in the regime of Hyper -Luminous IR-Galaxies (HyLIRGs, 
with Lfir ^ lO^^L©; e.g. 


Eisenhardt et al. 


20121. The corre¬ 
sponding star formation rate is estimated to be SFR ~ 3000 


M© yr (Sect. 4.2). HST/NICMOS imaging revealed tail¬ 


like features reminiscent of a merger system (Pentericci et al. 
|200I|| ). MRC 0152-209 is thus an ideal target to study the 
co-evolution of AGN and host galaxy. 


In [Emonts et al. ( |20IIa[ hereafter EMU) we pub¬ 
lished the detection of CO(I-O) emission in MRC 0152-209 
{L'qq ~ 6.6 X I0^° Kkm s~^pc^). We here present new 
high-resolution ATCA observations that allow us to map 
the spatial distribution of the CO(I-O). We will show that 
MRC 0152-209 is one of very few high-z galaxies with obser¬ 
vational evidence for tidal debris of cold molecular gas on a 
scale of ~60 kpc. Notably, the radio source is aligned with 
widespread CO(I-O) emission, suggesting a link between the 
cold molecular gas and the propagating radio jets. We also 
present new HST/WFPC2 imaging and WHT long-slit spec¬ 
troscopy to further investigate the violent merger episode 
that we witness in MRC 0152-209. 

Throughout this paper, we assume Ho=7I km s~^ 
Mpc“^, Dm=0.3 and r2A=0.7, correspondin g to an angul ar 
scale of 8.3 kpc arcsec”^ for MRC 0152-209 ( Wright|2006 ) ^ 


2 OBSERVATIONS 
2.1 CO(l-O) 

Observations of CO(I-O) were performed with the ATCA. 
We combined existing data taken on 25, 26 Aug. & 28, 29 
Sept 2010 with the most compact hybrid H75 and HI68 ar¬ 
ray configurations (baselines ranging from 31 — 192 meters; 
EMU) with new observations done on 23-25 June & 14- 
16 May 2011 using the inner 5 antennas of the 750B and 


^ See http://www.astro.ucla.edu/~wright/CosmoCalc.html 
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The Dragonfly Galaxy 3 


1.5B array configurations (covering baselines up to 1.3 km). 
Because phase stability decreases when baseline length in¬ 
creases or weather conditions deteriorate, we only used data 
taken under good weather conditions, i.e., data for which 
the root-mean-square path-length of the atmospheric seeing 
was an average [maximum] of 80 [200] /xm for the 750m and 
1.5k array configurations and an average [maximum] of 170 


[300] fim for the co mpact hybrid configurations (see Mid- 
delberg et al. 2006| lP] The total on-source integration time 


under these conditions was 40h (14.9h in 1.5k array, 12.6h 
in 750m array, 5.1h in H168 array and 7.8h in H75 array). 

The observational setup, calibration and data reduc¬ 
tion strategy was the same as described in [Ernonts et al.| 
( 2011a|b I. We used the Compact Array Broadband Back¬ 
end I Wilson et al.||20fT I, with two 2 GHz bands and 1 MHz 
channels. We centred both bands around zz = 39.46 GHz, 
which is the frequency of the redshifted CO(l-O) line (z/rest = 
115.2712 GHz). We took 2min scans on PKSB0130-171 to 
calibrate the phases and bandpass every 5-10 min, with the 
interval depending on the weather conditions and array con¬ 
figuration. Uranus was used to calibrate the absolute flux. 
Given that atmospheric phase instabilities increase signifi¬ 
cantly when going below 30° in elevation, most of our ob¬ 
servations were scheduled to avoid very low elevations. Be¬ 
cause the ATCA only has east-west baselines in its extended 
750m and 1.5k array configurations, this resulted in a lack of 
long-baseline uv-coverage at large hour angles, and thus an 
elongated beam with significantly worse spatial resolution 
in NS compared to EW direction. 

The data were calibrated in MIRIAD | Sault et al.|1995 1 
following EMU. We used a robustness parameter of 0.5 
(Briggs 19951 to image the data that are presented in this 
paper. This resulted in a beam-size of 4.0” x 1.3” (PA = 
-4.4°). The uncertainty in absolute flux calibration is ~30% 
(EMU). The final data set (covering ~ 15, 000 km s“^ with a 
full resolution of 7.6 km s“^ per 1 MHz channel) was binned 
to a channel width of 30 km s~^ and subsequently Hanning 
smoothed to an effective velocity resolution of 61 km s~^, 
with a noise level of cr = 0.146 mJy beam“^ chan“^. The 
data in this paper are presented in optical barycentric ve¬ 
locity relative to the redshift oi z = 1.9212 that we derived 
from the total GO(1-0) emission-line profile in EMU. 


2.2 HST imaging 

An optical image of MRG 0152-209 was taken with the 
Hubble Space Telescope (HST) Wide Field and Planetary 
Camera 2 (WFPC2) on 31 Oct 2000 (Project #8183, PI: 
G. Miley). The object was placed in the Planetary Gam- 
era frame, which has a pixel size of 0.0455 arcsec. The 
F555W broad-band filter (4790 - 6020 A) was used, which 
covers the near-UV rest-frame (1640 — 2062 A) for z = 1.92. 
The F555W filter also includes the redshifted HeHi 64 o and 
G HI] 1909 . However, from the spectrum that will be discussed 
in Sect. |3.3[ we estimate that the contribution of He H and 
GHI] to the total flux in the broad-band filter is <10%. 


^ Quoting |Middelberg et aL| ( |2006[ |: “The ATCA seeing monitor is 
a [2-dish] interferometer on a 230m east-west baseline that tracks 
the 30.48 GHz beacon of a geostationary communications satellite 
at elevation 60° ” 


We aligned the HST/WFPC2 image with the radio 
data by assuming that the central host galaxy is co-spatial 
with the peak of the GO (1-0) emission. This GO (1-0) peak- 
emission also coincides with the brightest component in the 


radio continuum from Pentericci et al.|2000 (Sect. [3^. The 
resulting HST astrometry agrees with the 3.6/im peak in 
existing Spitzer IRAG channel 1 data from [Seymour et al.| 
2007] (which has an estimated uncertainty in absolute as¬ 


trometry of ~0.3”). However, because the exact radio core 
is not known from existing radio continuum data, we con¬ 
servatively assume that the astrometric uncertainty in our 
HST overlay is ~0.5 arcsec. 

We also updated the astrometry of the near-IR 
HST/NICMOS F160W image from Pentericci et al. (20011 
by matching its nuclear features to the same features seen 
in the WFPC2 image. This NIGMOS image covers the 
4840 — 6230 A restframe emission at z = 1.9212. 


2.3 Spectroscopy 


A long-slit spectrum of MRC 0152-209 was taken on 25 July 
2012 using the Intermediate-dispersion Spectrograph and 
Imaging System (ISIS) on the William Herschel Telescope 
(WHT), La Palma. We used the R300B/R316R grating on 
the blue/red arm, the D5300 dichroic and the GG495 block¬ 
ing filter in the red arm to cut out second order blue light. 
This setup resulted in an effective wavelength coverage of 
3450 — 7270 A (1181 — 2490 A in the rest-frame), with a 
A-resolution of ~5A. The slit had a width of 1.3 arcsec. 
Given that the final CO(l-O) results were not available at 
the time of the WHT observations, the slit was positioned 
mid-way between the two 8.2 GHz radio continuum compo¬ 


nents imaged by Pentericci et al. (20001, at RA=01:54:55.74 


and dec=-20:40:26.42. Because of the low elevation at which 
the observations had to be done (airmass 2.4), we positioned 
the slit at the parallactic angle of PA = 320°. The on-source 
integration time was 3600 sec, divided into 6 x 600 sec con¬ 
secutive exposures. 

For the data reduction we used the Image Reduction 
and Analysis Facility (IRAF; [Tody 1993). A standard cal¬ 
ibration was done on each frame, including bias subtrac¬ 
tion, flatfielding and wavelength calibration. For the latter, 
we used an arc exposure taken at approximately the same 
time and elevation as the target spectra. The frames were 
then combined to remove cosmic rays and the background 
was subtracted. Four standard stars (Feige98, FeigellO, 
BD-l-253941 and Hiltner 102) were used for flux calibration. 

For the analysis of the spectra, we used the Starlink 
packages FIGARO and DIPSO. We binned consecutive rows 
of pixels to create a 1-D spectrum that contains the inte¬ 
grated flux across the full region were line-emission from 
MRC 0152-209 was detected. This spectrum was analysed by 
fitting Gaussian profiles to the emission lines. The spectrum 
is shown in the observed wavelength frame in this paper. 


3 RESULTS 

3.1 Stellar/HST morphology 

Figure shows the HST imaging of MRC 0152-209. We ar¬ 
gue that MRC 0152-209 is a system that is undergoing a 
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Figure 1. Left; HST/WFPC2 F555W image of MRC 0152-209. The four distinct features are described in the text. The inset on the 
top show a zoom-in of the central region #1. Right: HST/NICMOS F160W image at the same scale as the HST/WFPC2 image. The 
NICMOS data have been reproduced from Pentericci et al (2001), their Fig. 2. © AAS. Reproduced with permission. 


major merger, based on several distinct features: 1). A dou¬ 
ble nuclear component (separation 0.44" or 3.7 kpc). This 
double component is visible both in the WFPC2 and NIC¬ 
MOS imaging and could indicate the presence of either a 
close double nucleus or central dust feature; 2). A broad 
plume-like feature stretching 20 kpc (2.4") north of the cen¬ 
ter. From the WFPC2 image, this feature appears to be a 
tidal arm looping back onto the host galaxy; 3). A southern 
tidal arm/tail that is most prominent in the NICMOS image. 
In the middle of the tidal arm (at roughly 10 kpc distance 
from the center) there is an apparent NICMOS companion 
that is not visible in the WFPC2 image. Both this apparent 
companion and the middle part of the tidal arm therefore 
appear to be heavily obscured in the near-UV restframe; 4). 
A faint SW tail that ends in a region with enhanced emis¬ 
sion in both NICMOS and WFPC2, roughly 23 kpc (2.8") 
from the center. All these features are indicated in Fig.[2 


3.2 Molecular CO(l-O) gas 

In EMU we presented the detection of CO(l-O) emission 
in MRC0152-209, with Aco(i-o) ~ 6.6 x lO^^L©. Figures 
[^andj^show the distribution of CO(l-O) in MRC0152-209 
that we derive from our new high-resolution ATCA observa¬ 
tions. Part of the CO(l-O) emission appears to be associated 
with the bright optical emission of the central host galaxy 
(with ‘host galaxy’ defined in this paper as the region of 
bright optical emission that covers the double nuclear com¬ 
ponent in Fig.[^ i.e. region i^l). However, a large fraction of 
the CO(l-O) emission is clearly off-set from the host galaxy. 
This gas appears to coincide with several of the tidal fea¬ 
tures seen in the HST/WFPC2 and HST/NICMOS imaging, 
namely the inner parts of the looped northern arm (region 
#2) and the southern arm/tail (region #3). Following Pa- 
[padopoulos et al.|2008 (see also |Fomalont|199^ , we estimate 
that the astrometry of the CO emission-line peak in regions 
2 and 3 has an accuracy of 66rms ~ |(Obeam)(S/N)“^ ~ 
0.5"(NS) X 0.2"(EW), for ©beam ~ 4" x 1.3" and S/N - 4.5. 
This confirms that part of the CO(l-O) is indeed off-set by 


at least 20 kpc from the central host galaxy. There is also a 
weak tail of redshifted CO(l-O) that is not detected in the 
HST image and that stretches up to ~35 kpc (4") SE of the 
host galaxy (region “T” in Figs.[2]-[^. We argue that the 
CO(l-O) in regions 2, 3 and “T” is most likely tidal-debris 
of cold molecular gas that follows the tidal features visible in 
the HS'T imaging, and which is spread across a total extent of 
~60 kpc. Alternatively, some of the extended CO(l-O) com¬ 
ponents may identify gas-rich merging companion galaxies. 
We do not favor this scenario, because no identifiable candi¬ 
date galaxies are found at the location of the extended CO(l- 
0) in the HST/WFPC2 and NICMOS imaging. However, 
given the ~0.5" uncertainty in absolute astrometry between 
the ATCA and HST data, and the presence of a possible 
HST/NICMOS companion ~10kpc towards the South-East 
(Fig.0, we cannot completely rule out this possibility with 
our current data. 

Figure]^ shows the total intensity maps and spectra of 
the various CO(l-O) components. We derive CO luminosi¬ 
ties of the three main components (1, 2 and 3) from their 
emission-line spectra. We use this approach because the CO 
emission from these three components is not completely mu¬ 
tually independent (given that their spatial separation is on 
the order of one beam-size), and in the spectra any contri¬ 
bution from adjacent components can be more easily iden¬ 
tified than in the total intensity plots (Fig.|^- bottom). We 
estimate that Lqq = 2.7 ± 0.3 x 10^° K kms“^ pc^ is asso¬ 
ciated with the central host galaxy (the “blue” component 
#1 in Figji]), while L'co = 1-2 ± 0.3 and 1.4 ± 0.3 x 10^° 
K kms“^ pc^ are associated with the suggested tidal debris 
of the green (#2) and red (#3) component, respectively!^ 
Combined, these three components recover roughly 80% of 

^ Uncertainties quoted in this Section were calculated based on 
the FWZI of the Gaussian fits in Fig.^ (bottom), following Equa¬ 
tion 2 from [Emonts et al.| | |201^ based on |Sage||1990| |. They do 
not reflect the additional 30% error in absolute flux calibration 
(Sect. 113, because this is a scaling factor that affects all compo¬ 
nents in the same way. 
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l‘'54“56®0 55^8 55^6 55^4 55^8 55^6 55^4 55®8 55^6 55^4 55^8 55^6 55^4 

RA (J2000) RA (J2000) RA (J2000) RA (J2000) 


Figure 2. Channel maps of the CO(l-O) emission in MRC 0152-209. Contours of CO(l-O) emission are overlaid onto the HST/WFPC2 
image of Fig.[^ Contour levels: -3, -2 (dashed), 2, 3, 4, 5 (solid) xcr, with cr = 0.146 mjybeam“^ the root-mean-square (rms) noise level. 
The distinct regions indentified in the plot correspond to the same regions as in the HST image of Fig.[^ with “T” an additional tail-like 
feature that is only seen in CO (1-0). Note the highly elongated beam (ellipse) that we obtained by using the extended East-West array 
configurations of the ATCA (see Sect .|2.l| for details). 



Figure 3. Channel maps of the CO(l-O) emission, shown in RA vs. velocity. The data are binned to 1.2" in declination. Contour levels: 
-3, -2 (dashed), 2, 3, 4, 5 (solid) xcr, with a = 0.143 mjybeam“^ the rms noise level. The distinct regions indentified in the plot 
correspond to the same regions as in Fig.[^ 
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Figure 4. Top leR: Total intensity image of MRC0152-209 integrated over three distinct velocity ra^es (blue = -100—1-22 km s“^; 
green = 53—113 km s“^; red = 113 — 172 km s“^) and overlaid onto the HST/WFPC2 image of Fig.j^to within an accuracy of <0.5" 
(part of this figure is reproduced, with permission, from |Emonts et al.|2014] their Fig. 4). CO(l-O) contour levels are at 2.8, 3.5, 4.2, 4.9, 
5.6, 6.3(T level (with a = 12 mjy beam“^ X km s“^ for the blue component and 7.3 mjy beam“^ X km s“^ for the red/green components). 
The dashed line visualises the extrapolated radio axis. Right: Plots of the individual total intensity maps, as extracted across the three 
different velocity ranges. The 8.2 GHz radio continuum source from [Pentericci et al.| |2000| is overlaid in greyscale and thin contours 
(levels at 2, 7, 25 and 86% of the peak intensity). Bottom: CO(l-O) emission-line spectra in the four distinctive regions indicated in the 
total intensity plot on the top. The data are Hanning smoothed and the velocity is shown with respect to z = 1.9212 (Sect.j^. The solid 
blue, green and red lines represent Gaussian fits to the CO emission from regions 1, 2 and 3, respectively. From the flux covered by these 
Gaussians we derived ^CO for each region. The dashed colored lines indicate the estimated spill-over emission from adjacent regions that 
are separated by one beam-size or less (this emission is not included in the determination of L^q). The thin black line is the overall fit. 
[A color version of this Figure is available in the online edition of the paper.] 


the total CO(l-O) luminosity from the low-resolution obser¬ 
vations of EMU {Lqq ~ 6.6 ±0.8 x 10^° K kms“^ pc^; see 
also Emonts et al. ]m4l ). This suggests that another 20% of 
the CO emission {Lqq ~ 1.3 ± 0.9 x 10^° K kms“^ pc^) is 
associated with even more widespread gas, such as the green 
extension NW of component 2 and the tentative tail towards 
the SE. Therefore, ~40% of the overall CO(1-0) emission in 
MRC 0152-209 (i.e. L' = 2.7 ±0.3 x 10^° K kms“^ pc^) ap¬ 
pears to be associated with the central host galaxy, while 
~60% {L' = 3.9 ± 1.0 X 10^° K kms“^ pc^) is found across 
the widespread tidal debris. Fig.also shows that the tidal 
CO(l-O) is found at higher velocities than the CO(l-O) in 
the central host galaxy. 


The 8.2 GHz radio continuum source appears to be 
aligned with the CO-emitting gas (Fig.[^- right). In par¬ 
ticular the alignment with the widespread CO(l-O) emission 
NW of the host galaxy is striking. Pentericci et al. 120001 


classify the radio source as a small double-lobed radio source 
with two steep-spectrum hot-spots (af ^ = —1.4 for the NW 


and -1.7 for the SE hot-spot), which cover a total linear ex¬ 
tent of ~18 kpc (2.2"). We are currently in the process of 
imaging the structure of the radio source at much higher res¬ 
olution with new observations. If confirmed that this is a jet, 
the alignment with the widespread CO(1-0) emission NW of 
the nucleus would suggests that there is a direct physical link 
between the cold molecular gas and the propagating radio 
source. We will discuss this further in Sect. 14.31 


3.3 Near-UV emission-line gas 

Figure shows the WHT optical long-slit spectrum of 
MRC 0152-209, which covers 1181 — 2490 A in the near-UV 
rest-frame. We clearly detect the emission lines of Lya, CIV, 
Hell and NelV, with a tentative detection of CH. Our mea¬ 
sured Lya flux is in good agreement with the value quoted 
by McCarthy et al.|(1991| . The CHI] line, as reported by 
McCarthy et al.j(1991|, fell on a prominent broad sky-line 


residual in our data (which we blanked-out in Fig.j^, hence 
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Figure 5. WHT long-slit spectrum of MRC0152-209. The bottom-left plot visualizes the slit alignment (see Sect. |2.3| for details). The 
lines of Lya, CIV, Hell and NelV are clearly detected, with a tentative detection of CII. The red line in the insets shows a Gaussian 
fit to the line-profiles. 


Table 1. Optical spectroscopy. Observed redshift z, FWHM (corrected for instrumental broadening), line flux/luminosity and relative 
strength with respect to Lya for the various emission lines detected in the optical spectrum. For completion, the z and FWHM of the 
CO(l-O) line emission form EMU, based on the integrated emission-line profile, are also included. 


Line 

z 

FWHM (km s"!) 

Siine (10 l®ergs Wm 2) 

Liine (lO^^^ergs i) 

R (Suya/Sline) 

Lva 

1.9215 ±0.0007 

1301 ±187 

30.5 ±2.8 

81.8 ±7.5 

1.0 

1215.7A 

CIV 

1.922 ±0.02 

900 ± 200 

2.56 ±0.24 

6.87 ±0.64 

11.9 

Hell 

1.9214 ±0.0007 

687 ±147 

1.23 ±0.20 

3.30 ±0.54 

24.8 

1640.4A 

NelV 

1.9211 ±0.0007 

839 ±148 

0.88 ±0.14 

2.36 ±0.38 

34.7 

2424.4A 

CO(l-0)ii5.3 GHz 

1.9212 ±0.0002 

397 ±37 

- 

- 

- 


could not be detected. The properties of these emission lines 
are summarised in Tabled 

The emission lines are all well approximated by a sin¬ 
gle Gaussian profile (Fig.[^. The FWHM of the emission 
lines (corrected for instrumental broadening) varies from 
1300 km s“^ for the strong Lya line to ~700 km s ^ for 
Hell, compared to FWHM~400 km s“^ for the integrated 
CO(l-O) profile (EMll)[^The relatively broad Lya profile 
can partially be explained by resonant scattering. Never¬ 
theless, even for the other emission lines of ionised gas, the 
line-widths are roughly twice that of the CO(l-O) emission. 
This is similar to what has been found when comparing the 
[O hi] and CO kinematics in low-z QSOs, and has been as¬ 
cribed to a different geometry and spatial extent for the 


^ The CO(l-O) profile from EMU was extracted across a larger 
region than the WHT spectrum. Fig.[^ shows that the near-UV 
emission from the central host galaxy was included in the slit, 
but that some of the extended emission may have been missed. 
Therefore, the FWHM of the near-UV lines are firm lower limits. 


ionised and molecular gas phase, as well as a higher sensi¬ 
tivity of the ionised gas to non-gravitational motions, such 
as outflows ( Shields et al.|2006 Villar-Martfn et al.|2013b |. 
No clear absorption features are observed in any of the emis¬ 
sion lines. We note, however, that the Lya line falls at the 
edge of our spectrum, where a significant drop in sensitivity 
resulted in an increased noise level. 

To within the uncertainty, the redshift derived from 
each of the optical emission lines is in excellent agreement 
with z = 1.9212 ± 0.0002 derived from the integrated CO(l- 
0) profile in EMU. 


4 DISCUSSION: THE ‘DRAGONFLY GALAXY’ 

We presented high-resolution CO(1-0) data, complemented 
with HST/WFPC2 imaging and WHT spectroscopy, of the 
merger system MRC 152-209 (z = 1.92). Based on the opti¬ 
cal morphology of the galaxy in our HST/WFPC2 image 
(Fig.0, we will here adopt the name ‘Dragonfly Galaxy’. 
In this Section we will further discuss the CO(1-0) proper- 
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ties (Sect. E3, compare those to the extreme star formation 
rates (Sect. HI and investigate AGN feedback (Sect. [43| . 
This will further reveal the evolutionary stage of this ob¬ 
ject. 


dio galaxy. This is in agreement with the claim by |Ivison| 


et al. (2012) that violent gas-rich galaxy interactions are 


ubiquitous among starbursting HzRGs. 


4.1 Molecular gas & fueling 

We saw in Sect. |3.2| that a large amount of cold molecular 
gas (Too 3.9 ± 1.0 x 10^° Kkm s“^ pc^) appears to re¬ 
side in large-scale tidal features that stretch several tens 
of kpc from the center of the Dragonfly Galaxy. We as¬ 
sume a standard conversion factor of Mh 2 /T'qo = 0.8 M© 
(K kms“^ pc^)“^ to translate our CO(l-O) estimates into 
molecular gas masses. This agrees with Mh2/T'co found in 
ultra-luminous infra-red galaxies ( Downes &: Solomon] 1998' I, 
as well as high-z sub-millimeter and starforming galaxies 


4.2 Star formation & evolution 


'acconi et al. 2008 Stark et al. 20081, although we dis¬ 


cussed in EMU that this conversion has a large uncertainty. 
This means that Mh 2 ~ 3.1 ± 0.8 x 10^° Mq is associated 
with the widespread gas reservoir, most likely consisting of 
tidal debris, that surrounds the Dragonfly Galaxy. An ad¬ 
ditional Mh 2 ~ 2.2 ± 0.2 X 10^° M© is found in the central 
host galaxy. 

The GO(l-O) in the tidal debris appears to be spread 
across a total extent of roughly 60 kpc. This is larger than 
the extent of tidal CO features observed in other high- 
z galaxies. For example, observations with the Plateau de 
Bure Interferometer showed CO(3-2) emission on scales of 


8 — 20 kpc in two high-z SMGs (Tacconi et al. 20081, while 


CO(2-l) observations with the Very Large Array showed 


tidal debris across 5 kpc in a z = 4.4 QSO (Riechers et al. 


20081. In fact, detecting cold gas on scales of even a few 


tens of kpc is still rare at high-z. While ALMA observations 
of the bright atomic [Clljiss^im have not detected any ex¬ 
tended cold gas (e.g.|Wang et al. 2013 [De Breuck et al.|2014 


[Riechers et al.|2014 1, jCicone et (20151 recently reported 

[CII] emission spread out to ~30 kpc from the center in a 
z = 6.4 QSO. The latter detection was obtained with the 
most compact conhgurations of the Plateau de Bure Inter¬ 
ferometer. We therefore argue that sensitive short-baseline 
observations are key to detecting such very extended gas 
reservoirs. In this respect, the relatively large number of 
short (<200m) baselines that were included in our ATCA 
observations likely provided the sensitivity for tracing CO(l- 
0) emission on tens of kpc scales. 

On the other hand, the very extended molecular gas 
reservoir around the Dragonfly Galaxy may also serve as 
an extreme example that gas-rich galaxy mergers and in¬ 
teractions may be a common triggering mechanism for the 
powerful radio sources in HzRGs. This would be analogue to 
what is believed to occur for only the most powerful (high- 
excitation) radio sources at low- and intermediate-z, which 
are thought to be fueled by cold gas deposited after a ma¬ 
jor merger event (e.g. Hardcastle et al. ][20^ [Emonts et al. 


20I0II [Smolcic fc Riechers[[20II| [Janssen et al.[[20I2f[R 


Almeida et al. 2012 Mao et ah] 20141. Of course, our cur¬ 

rent work is based on a single object and cannot exclude 
that other possible mechanisms may provide fuel for star- 


formation and AGN activity in HzRGs (see e.g. Hatch et al. 
20131. Nevertheless, our CO(l-O) results on the Dragonfly 


Galaxy provide direct observational evidence of a gas-rich 
(‘wet’) major-merger associated with a powerful high-z ra- 


Drouart et al. (2014) show that the Dragonfly Galaxy is 


the most IR-luminous HzRG known in the southern hemi¬ 
sphere, with Lir( 8 - 1000 Mm) ~ 2.2 x IO^^Lq, which is in the 
regime of HyLIRGs. About 80% of the IR emission is es¬ 
timated to arise from the dust-enshrouded starburst com¬ 
ponent (i.e. Lfir-sb ~ 1-8 x lO^^L©), while the remain¬ 
ing 20% comes from the buried AGN. Following [Drouart] 
et al. ( [2014 ), this corresponds to a star formation rate of 
SFR~3000 MQyr“^. The SFR is a factor of a few higher 
than any other HzRG at z < 2.5 studied by [Drouart et al~| 
(2014), and only rivaled by three known northern HzRGs at 


z > 3.5. For example, the estimated starburst-IR luminosity 
of the Dragonfly Galaxy is more than 3x that of the mas¬ 
sive Spiderweb Galaxy at z = 2.1 and comparable to that 
of 4C41.17 at z = 3.6 ( [Drouart et al.[[2014 |. The Dragonfly 
Galaxy is also among the most massive HzRGs known at 
any z, with M* =5.8 x 10^^ M© ( De Breuck et al^|2010| . 

As shown in Fig.j^ and Table strong Lya emission 
is detected in the Dragonfly Galaxy. The total Lya lumi¬ 
nosity covered by the 1.3” wide slit is Lnya = 8.2 ± 0.8 x 
10"^^ erg s“^. This is within the range of Lcya that Rottgering 


et al. (19971 detect among a sample of 30 HzRGs at z > 1.9. 


Evidence that part of the Lcya in the Dragonfly Galaxy orig¬ 
inates from enhanced star formation comes from the high 
Lya/Hell emission-line ratio of 25 (Table [^. According to 
[Villar-Martfn et al. (2007), a line ratio of Lya/He 11^15 
indicates that simple AGN photo-ionisation models are no 
longer able to explain the Lya content, and instead a sig¬ 
nificant fraction of the Lya emission is ionised by enhanced 
star formationj^The Lya/Hell = 25 found in the Dragonfly 
Galaxy is only matched by three other HzRGs out of a sam¬ 
ple of 54 (4G41.17, PKSB124 3-f036 and TN J0205-f22 42), 
which are all three at z > 3.5 (Villar-Martfn et al.j2007|. In 


addition, the corresponding LLya/LiR-SB ~ 1-2 x 10" “ is an 
order of magnitude below the estimated LLya/LpiR ~ 0.03 
that Humphrey et al. ( 2011| derive from the expected in¬ 
trinsic Lya luminosity of massive starbursts. Gombined with 
the fact that Lya has a contribution from gas that is pho- 
toionised by the AGN, it is clear that the Lya emission 
due to stellar ionisation is much fainter than expected for 
a system with such intense star formation. This means that 
the starburst in the Dragonfly Galaxy must be very heav¬ 
ily obscured. These Lya results support that the Dragonfly 
Galaxy is an extreme z ~ 2 starbursting system, in which 
the enhanced star formation is presumably triggered by the 
merger event. 

When comparing the starburst IR luminosity to the CO 


^ An enhancement in the collisionally excited Lya and CIV lines 
with respect to He II can also occur from high gas densities (>1000 
cm“^), but this would imply that the Lya emission is mostly un¬ 
obscured, which is not likely for the Dragonfly Galaxy. To confirm 
that enhanced star formation causes the large Lya/Hell ratio, 
deeper optical spectra are required to search for stellar absorp¬ 
tion features similar to those found in 4C 41.17 ([Dey et al.|19Tf|. 
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The Dragonfly Galaxy 9 


luminosity, Lpir/Lqo ~ 300 Lq (K km/s for the 

Dragonfly Galaxy. This is comparable to what is observed in 
the IR-brightest high-z sub-millimetre galaxies (SMGs; e.g. 


Greve et al.|2005]|Genzel et al.| 20 T 0 l|Ivison et al.|2011[|Both^ 


well et al.|2013 1. The specific star formation rate of the Drag¬ 


onfly Galaxy, sSFR = SFR/M* =5.2 Gyr ^, is also compara¬ 
ble to the sSFR~ 1 — 10 Gyr“^ estimated for SMGs and high- 


z dust-obscured galaxies (DOGs) by Bussmann et al. (2012 1 . 
These estimates are higher than those obtained by cosmolog¬ 
ical simulations of SMGs undergoing smooth gas accretion 


(e.g. Dave et al. 20101, leading Bussmann et al. (20121 to 
conclude that the high sSFRs in these systems are driven by 
galaxy mergers. For the Dragonfly Galaxy, with SFR~ 3000 
Mq yr“^ and Mh 2 ~ 5.3 x 10^°Mq, the estimated mini¬ 
mum gas-depletion time-scale is tdepi = Mh 2 /SFR ~ 18 
Myr (assuming that the complete reservoir of cold molecu¬ 
lar gas is efficiently converted into new stars). If the bulk 
of the star formation occurs in the central host galaxy, with 
Mh 2 ~ 2.2 X 10^° M© of cold molecular gas available, this 
would be tdepi ~ 7 Myr. These estimates are at the low end 
of the gas-depletion time-scales typically derived for high-z 


SMGs ([Greve et al.||2005| Tacconi et al. [2008 Ivison et al. 


2011 


[Riechers et al.||2011c[ but see also [Swinbank et ~ 


2006 1 . They are also about an order of magnitude shorter 


than the gas depletion time-scales derived for a sample of 
starforming galaxies (SFGs) at z = 1 — 3 (Genzel et al. 2010| , 
or for a similar sample of IR-selected SFGs galaxies at z ~ 2, 
which dominate the ULIRG-mode at high-z through gas- 


accretion processes other than major mergers (Daddi et al. 
2007 Rodighiero et al)]|2011 l. This IS again consistent with 
the interpretation that the current HyLIRG-phase of the 
Dragonfly Galaxy is a short period of maximum starburst 
activity triggered by a major gas-rich merger. As such, the 
Dragonfly Galaxy resembles the well-studied gas-rich galaxy 
group BRI 1202-0725 at z = 4.7, which contains a merg¬ 
ing dust-obscured SMG and FIR-hyperluminous quasar (e.g. 
Wagg et al. 2012| Salome et al. 2012| Garilli et al. 20131. 
However, with a separation of ~26 kpc between the SMG 
and quasar in BRI 1202-0725, the notable difference is that 
the Dragonfly Galaxy appears in a much more advanced 
stage of the merger. 

Interestingly, the minimum time-scale for the star for¬ 
mation to deplete the cold gas is of the same order as the 


expected lifetime of the radio source (e.g. Liu et al. 1992 
Blundell fc Rawlings||200^. This short depletion time-scale 


implies that —in principle— the star formation itself could 
be the dominant mechanism for clearing the central host 
galaxy of cold, self-gravitating molecular gas, without the 
explicit need for strong negative feedback from the quasar 


2005 


nucleus (as often invoked for high-z QSOs; e.g. Springel et al. 
Di Matteo et al.]|2005 1. This short and vigorous star- 


burst episode may perhaps indicate that we are witnessing 
a crucial phase of transformation that will quickly turn the 
Dragonfly Galaxy into a massive high-z quiescent elliptical 


(e.g. van de Sande et al. 20131, once most of the cold gas 


has been consumed. 

Still, despite the estimated short tdepi for the current 
starburst event (which is most likely concentrated towards 
the central host galaxy), it is likely that the star formation 
rates within the widespread tidal debris are substantially 
lower than 3000 M 0 yr“^. Given that roughly 60% of the 
cold molecular gas is spread on scales of this tidal debris. 
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it is likely that part of this gas reservoir will remain intact 
for much longer than tdepi and may trigger future, perhaps 
less intense, episodes of star formation, once the cold gas is 
re-accreted back onto the host galaxy. 


4.3 Radio source & feedback 


The Dragonfly Galaxy has a fairly small radio source 
(2.2"/18kpc in diameter; Pentericci et al.||2000 |. On scales 
of a few tens of kpc, there is an apparent alignment between 
part of the widespread GO(I-O) emission and the radio-jet 
axis (Fig.|^. This is most notable for the GO-emission NW 
of the host galaxy. Unless this is merely a chance align¬ 
ment, it suggests that there is a physical link between the 
propagating radio source and the widespread molecular gas. 
Intriguing support for this scenario comes from the recent 
discovery of cold molecular CO(l-O) reservoirs found in the 
halo environment of several HzRGs, with the GO(l-O) emis¬ 
sion aligned along the radio axis and found just outside the 
brightest edge of the radio source ( Emonts et al.|2014 ). 

The widespread CO (1-0) emission in the Dragonfly 
Galaxy suggests that it may be common for cold gas and 
dust to be expelled as tidal debris during the early evolu¬ 
tion of high-z radio galaxies. As the jets expand into this 
cold gas and dust, on the one hand their radio flux may be 
boosted, which may potentially introduce a bias in identify¬ 
ing these systems as HzRGs (jBarthel Sz Arnaud|1996 Tad- 


hunter et al. ]MI| . On the other hand, it is also likely that 


slow shocks will develop in the multi-phase medium, leading 
to compression of the atomic and molecular phases of the 
gas ( Guillard et al.|[^009[ ). Small amounts of dust surviving 
in the shocked gas may catalyze the formation of molecules. 
The cooling provided by (predominately) H 2 may lead to 
large quantities of cold molecular gas (Guillard et al.[20l0| 


as observed along the radio jets of several HzRGs ( Klamer 
et al. 2004| Emonts et al. 20141. Thus, the tidal stripping 


of gas and dust during a galaxy merger may possibly play 
a key role in enabling the formation of cold molecular gas 
along the radio jet as it expands. 

Potentially, this scenario may even trigger jet-induced 
star formation, which has been invoked as a possible ex¬ 
planation for alignments seen between radio sources and 
UV/optical emission in HzRGs (McCarthy et al. 1987| 


Chambers et al. TT9871 see also follow-up work by e.g. |Stein-| 
bring||2014| Zinn et al.||2013 Karouzos et al.|[20I4 l. In this 


respect, it is interested to note that the radio source in the 
Dragonfly Galaxy also aligns with the main axis of the host 
galaxy in the HST imaging of Fig.[2 


5 CONCLUSIONS 

We have provided evidence that the Dragonfly Galaxy 
(MRC 0152-209) is a high-z ‘wet’ merger, observed roughly 
3.5 Gyr after the Big Bang. With its extreme star formation 
and evolving radio source, the Dragonfly Galaxy allows us 
to witness a short but crucial stage in the early evolution 
of massive galaxies, during which the mass build-up is ex¬ 
pected to be at a maximum for a period of a few tens of 
Myr. This paper highlighted the crucial role that the cold 
molecular gas is expected to play in this. 

The main results of this study are: 
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i). HST imaging shows the presence of large-scale tidal fea¬ 
tures that identify the Dragonfly Galaxy as a merger. 

a). ATCA CO(l-O) mapping reveals that cold molecular gas 
is associated with the tidal debris, classifying the Dragonfly 
Galaxy as a ‘wet’ merger. We conclude that: 

- the cold molecular gas is distributed across a total extent 
of ~60 kpc, making the Dragonfly Galaxy one of very few 
high-z objects in which the CO has been imaged on these 
scales; 

- roughly 60% of the CO(l-O) emission (Mh 2 ~ 3x 10^° M©) 
appears to be associated with this widespread tidal debris, 
while the remaining 40% (Mh 2 ' ^ 2 X 10^° M©) is found 
in the central host galaxy; 

- the minimum gas depletion time-scale is a short 7—18 
Myr (although the gas reservoir at large scales may trigger 
future star formation episodes, once re-accreted back onto 
the galaxy); 

- the Dragonfly Galaxy may be a high-z analogue to the 
most powerful radio sources at low-z, whose activity is 
believed to be fueled by cold gas deposited after a galaxy 
merger. 

Hi). Consistent with the very high star formation rates 
(SFR~3000 M© yr“^) found in previous IR studies of 
the Dragonfly galaxy, our new WHT spectrum confirm 
the presence of enhanced and heavily obscured star forma¬ 
tion through a high Lya/Hell emission-line ratio and low 

Lhya/LlR-SB- 

iv). The small (18 kpc) radio source appears to be aligned 
with widespread CO (1-0) emission, suggesting a physical 
link between the propagating radio jets and the cold gas 
reservoir. We propose a scenario where tidally ejected gas 
and dust that is overrun by a radio jet may undergo com¬ 
pression by slow shocks, resulting in the formation and cool¬ 
ing of molecular gas along the radio axis. 

Concluding, as the most IR-luminous southern HzRG, the 
Dragonfly Galaxy provides us with a unique opportunity to 
further investigate the most active stages of galaxy evolution 
in the Early Universe. 
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